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Abstract 

The properties of a generalized version of the Borel Transform in 
infrared unstable theories with dynamical mass generation are studied. 
The reconstruction of the nonperturbative structure is unambiguous 
in this version. Various methods for extracting the singularity struc- 
ture of the Borel Transform for lattice formulations of such theories are 
explored, and illustrated explicitly with the O(N) sigma model. The 
status of the first infrared renormalon in QCD is discussed. The feasi- 
bility of a proposed technique for analytically continuing from the left 
hand Borel plane (where nonperturbative information is available via 
simulation of lattice field theory) to the positive real axis is examined 
using the sigma model. 
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1 INTRODUCTION 



The problem of divergence of (renormalized) perturbation theory, which can 
be traced back to the seminal work of Dyson in quantum electrodynamics 
[Q], is now recognized to lie at the core of any precise analytic understanding 
of the nonperturbative structure of field theory. The need to face this prob- 
lem squarely is particularly apparent in massless quantum chromodynamics 
(QCD), where the underlying dynamics is specified by a single parameter N 
(for gauge group SU(N)). Dimensionless ratios of hadron masses in this the- 
ory depend only on N, which is of order unity (N =3) in the physical world. 
Apart from the technically intractable 1/N expansion, there is no natural in- 
trinsic expansion parameter for the spectral and other low-energy properties 
of the theory. If one insists on computing such quantities in a conventional 
weak-coupling expansion, they either vanish formally to all orders, or yield 
asymptotic expansions in terms of a running coupling which is naturally of 
order unity. Such expansions are quantitatively useless in the absence of a 
reliable "resummation" procedure. 

Certain superrenormalizable field-theories (e.g <f>\ 2 3 in the unbroken phase) 
have been shown rigorously to possess Borel-summable perturbative 

expansions, so that the full content of the theory (at least, information- 
theoretically) is exhausted by perturbation theory. The spontaneously bro- 
ken phases of these theories is typically not Borel-summable, but even here 
an optimized reorganization of the perturbation theory can be shown rigor- 
ously [|J to converge to the exact partition function, provide the theory is 
formulated at finite volume. The failure of Borel-summability is far more 
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severe in QCD ||, and no systematic analytic procedure is known, even in 
principle, whereby the full nonperturbative structure of the theory could be 
obtained on the basis of purely perturbative information. 

The divergence of the perturbation theory in QCD is far more than a 
merely technical embarassment, as has recently been emphasized by Mueller 
|J. It fundamentally limits our ability to reliably compute important correc- 
tions (such as higher twist effects) to the vast phenomenology of high energy 
processes described by perturbative QCD. Specifically, in those cases where 
one attempts to weld analytic perturbative with numerical or phenomenologi- 
cal nonperturbative estimates for the same process, the precision of the result 
is clouded by unavoidable resummation ambiguities on the perturbative side. 
Consider, for example, a renormalization group controlled amplitude H{Q 2 ) 
in QCD, expressible as a formal (divergent asymptotic) series in the running 
coupling a(Q 2 ): 

n(Q 2 )^Y,c n a n (Q 2 ) (1) 

n 

Quite generically, one may identify subsets of graphs in perturbation theory 
which contribute factorial growth to the dimensionless coefficients c n at large 
n 

c n ~ (-)Vn! (2) 
V 

where bo, b\ = hj ^- are the first two coefficients of the beta function, and p 
is a positive integer. Using Stirling's approximation, the ra-th term in the 
asymptotic series (0) may be written 

n!(^) n n 7 a(Q 2 )" ~ e »M»»Wp)-r» ( 3 ) 
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which is minimal at n m i n ~ feoa (Q2) — P 

lnfr. The error in an asymptotic 
expansion is typically at least as large (it may of course be much larger !) as 
the smallest term, here 



which is power suppressed at large Q in just the way we expect for a higher 
twist term. The moral is clear: higher twist effects are not precisely calculable 
in the absence of a reliable resummation procedure. The same arguments 
in fact apply in a variety of other important situations: for heavy quark 
expansions, where 1/M P effects (M=heavy quark mass) may be confused 
with resummation ambiguities arising in a static-quark quantity evaluated 
via the renormalization group in terms of a perturbative expansion in a(M), 
or in estimating analytically power corrections to lattice quantities in the 
lattice spacing a. 

The conventional Borel transform approach (for a thorough review, see 
[0]) attempts to reconstruct the full nonperturbative structure of the the- 
ory from knowledge of the behavior of weak-coupling perturbation theory 
at large orders. There are two serious problems with such an approach in 
infrared unstable theories like QCD. First, it is becoming increasingly ap- 
parent that the large-order behavior is exceedingly complex ||. It is very 
difficult, even restricting oneself to a strictly limited subset of graphs, to 
precisely identify the terms which give the truly dominant behavior at large 
order, due to the subtle interplay of combinatoric and kinematic effects. Sec- 
ondly, in a non-Borel-summable theory (and any infrared unstable theory 
with dynamical mass generation falls into this category) there is in any case 




(4) 
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no reliable reconstruction theorem based on the usual perturbative Borel 
transform, which develops singularities (renormalons) on the positive real 
axis of the Borel variable s, i.e on the Borel integral contour. Fortunately, 
an alternative formulation of the Borel transform exists, first discussed in 
QCD by t'Hooft [|, and later by Crutchfield @ and David 0. This form 



of the transform is closely connected to the "naive" perturbative definition 
in Borel-summable theories, but can be used to give a precise reconstruction 
theorem in a much wider class of theories. The singularity structure of this 
transform is unambiguously related to the exact nonperturbative structure of 
the full theory. But by the same token, this structure is not directly accessi- 
ble in conventional perturbation theory, but must be studied in an explicitly 
nonperturbative way. 

The primary objective of this paper is to suggest that useful informa- 
tion about the renormalon structure (i.e positive s singularities) in infrared 
unstable theories can indeed be obtained using nonperturbative techniques 
(such as the large N expansion, or lattice theory). In Section 2, we review 
very quickly the conventional perturbative Borel transform B(s), and ex- 
plain its relation to the generalized transform B(s) and the density of states 
function B(s) = disoB(s). A classification of renormalon singularities into 
those which are already present at finite volume ("Type 1") and those which 
are only strictly speaking present at infinite volume ("Type 2") is given. 
The basic interconnections are further illustrated using spin models in 1 and 
2 dimensions (the nonlinear sigma model) in Section 3. In Section 4 we 
explain the perturbative camouflage of renormalon singularity structure in 
disconnected quantities such as the partition function (vacuum- vacuum am- 
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plitude). The status of the first infrared renormalon singularity in QCD 
(recently raised by Brown, Yaffe and Zhai [11]; see also |12|]) is discussed in 
Section 4. It is shown that certain renormalon singularities of the full Borel 
function B(s) are immune to infrared cutoffs in ordinary perturbation theory 
which remove the usual IR renormalons. Moreover, the presence of such a 
singularity appears to be a generic feature once dynamical mass generation 
is assumed. In Section 5 we study the sigma model at large N on the lattice 
in order to establish the numerical feasibility of extracting renormalon sin- 
gularities by analytic continuation of nonperturbatively computed (say, by 
Monte Carlo simulation) Euclidean amplitudes. Section 6 summarizes our 
conclusions and indicates directions for further study. 



2 Borel Transform Technology 

The archetypal example of a Borel-summable perturbation theory is the non- 
Gaussian integral 

Z(f) = f +0 ° e-7 v(x) dx, V(x)=x 2 + x 4 (5) 

= / e~ x ~ fx dx (6) 
Jo 

- E(-/r r(2r ?^ E c ^ w+1 ( ? ) 

n=0 n - 

Evidently the coefficients c n in the asymptotic expansion of Z(f) for / — > + 
have the large order behavior c n ~ (— l) n A n n\ . Define the "naive" Borel 
transform 

B(s) = E ^ (8) 
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so that Z (f) = J +o ° e~ s lf B(s)ds and Z(f) have the same asymptotic ex- 
pansion for small /. In this case the analytic properties of Z(f) are such 
(see Ref |7j for precise conditions) that Z = Z and knowledge of the "naive" 
transform leads to a precise reconstruction of Z for all /. We shall actually 
be using a more powerful version of the Borel method in most of this pa- 
per, which will circumvent the need for theorems of this type whcih certify 
the preceding procedure. The shortcomings of B(s) for non-Borel-summable 
theories, where the sign oscillation in c n is absent in some uncancelled set 
of contributions at large order (not necessarily the dominant onesl) are im- 
mediately apparent: B(s) develops singularities on the positive real s axis, 
rendering the reconstruction integral ambiguous (if the singularities are non- 
integrable). The most straightforward way in which such singularities can 
arise is illustrated in Fig. [1], where the "action" V(x) possesses a secondary 
nonperturbative extremum. 

To expose the singularity in B(s) write [13 

Z(f) = J e- v(x)/f dx (9) 

z'+OO r 

= / dse~ s/f / dx5(s - V{x)) (10) 

f+oo _ 

= / dse- s/f B(s) (11) 
Jo 

where 

B(s) = J dx6(s - V(x)) = E — 1^ (12) 

V(xi)=s I v */ I 

If V has a nontrivial extremum where V'(xi) — 0, Sj = V(xi) > 0, we will 
have necessarily B(si) = oo, a singularity (in this case, an integrable square- 
root singularity, as B ~ 1/ yjs — Sj) for positive s. To the extent that the 
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Figure 1: NonGaussian Action leading to singular B(s) 

interchange of integrations carried out in Eq@ is legal, the reconstruction 
of Z from B must be precise. In fact B is intimately related to the more 
general Borel function to be introduced below. In general, however, B ^ B\ 
Of course, the asymptotic expansions for small s of B, B, which determine 
the asymptotic expansion of Z(f) for small /, must agree. As pointed out 
originally by t'Hooft []5J, the existence of Euclidean extrema of the action 
(instantons) in QCD implies in a similar way the presence of singularities 
in B(s) at s = S inst . Moreover, the location of the singularities (though 
not their strength) is universal (see (0): once a singularity appears in any 
amplitude of the theory, it is expected to propagate to all others. 

For a class of non-Borel-summable theories (the archetype is the dou- 
ble well anharmonic oscillator) an optimized reorganization of perturbation 
theory can be shown to converge with exponential rapidity Q to the exact 
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partition function. These convergence proofs apply only when the theory is 
formulated at finite volume, however, limiting the practical usefulness of such 
methods in higher dimensional field theories. The approach taken here will 
be to employ a generalized Borel Transform which allows a "bulletproof re- 
construction of the full theory, avoids resummation ambiguities, and is more 
closely tied to the path integral formulation of the theory |§. Our ultimate 
objective is to explore the potential of the presently available nonperturbative 
lattice techniques for computing the singularity structure of this transform. 

Consider a typical infrared-unstable (asymptotically free) theory, such as 
the O(N) sigma model in 2 dimensions, or QCD in 4 dimensions, defined by 
a functional integral 



where is a generic set of fields, S [(/>]) the Euclidean action, and / a bare 
coupling constant. A compact lattice formulation (/ T>(f) = 1) of the theory 
is supposed, so that a UV cutoff is always present, and we have logarithmic 
asymptotic freedom 



where m ~ (physical correlation length) -2 , A ~ (short distance cutoff) -2 . 
For both sigma models and QCD on a finite lattice, < S[<f>] < S mauX is a 
bounded continuous function of the fields, and the integral / Vcp is compact. 
Define the generalized Borel Transform B(s) as the Laplace transform with 
respect to inverse coupling, as follows 




(13) 



/ 



1 



A >> m 



(14) 





(15) 
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The integral is convergent at large / for all s (as e s ^ Z(f) — > 1) and for 
Re(s)< as / — > + . Thus B(s) is left-half-plane analytic, with a cut on the 
positive real axis for < s < S ma _ x . Indeed, inserting the functional integral 
expression (|T3| ) in flI5|), and defining x = j, 

B(s) = jVtpj e x{s - s[rp]) dx (16) 

= l v *mrs (17) 

which is manifestly analytic in the s-plane cut along the positive real axis 
between s = and s = S mSuX . The discontinuity of B(s) across this cut is 

lmB(s = SR + te) = nJ V<j)6(s R - S[(f>]) = nB(s R ) (18) 

The physical interpretation of this discontinuity is clear: B(s) is the density 
of configurations in the functional integral with action equal to s (cf. Eqfll^)). 
Singularities of B are expected wherever there is a sharp local enhancement 
in the number of configurations, for example wherever the action has a local 
extremum. 

The rigorous reconstruction of Z(f) from B is trivial in this framework. 
The vertical contour of the inverse Laplace transform may be wrapped around 
the cut on the positive s axis to give an integral over the discontinuity of 
B(s), i.e. over B. More directly, we may simply observe that it follows from 
EqQID that 

Z(f) = / B(s)e- s/f ds (19) 



o 



This reconstruction is exact: there are no routing ambiguities in the s- 
integral. Of course, B is in general a distribution (containing potentially S 
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singularities), not obtainable purely from knowledge of weak-coupling pertur- 
bation theory. Instead, it must be computed nonperturbatively by analytic 
continuation of B(s), computed nonperturbatively in the region Re(s) < 0, 
where the defining integral ( |T5f ) exists. It should be noted at this point 
that B(s), in contrast to B(s) may have singularities in the left-half-plane 
(corresponding to non-principal-sheet singularities of B(s)): the ultraviolet 
renormalons of the renormalized perturbation theory in QCD are of this type. 

It is convenient to distinguish between two types of singularities appear- 
ing in B(s). Exact local extrema of the Euclidean action may exist even 
for systems with a finite number of degrees of freedom (e.g. field theories 
formulated on a finite lattice). Such extrema give rise to actual singularities 
(henceforth called "Type 1" ) of B even at finite volume, when the theory is 
cutoff in the infrared. For example, one can find lattice analogs of the instan- 
ton solutions of the classical 2D 0(3) sigma model using a simple numerical 
annealing procedure to solve the lattice Euclidean field equations. Other 
singularities of B only appear at infinite volume, and will be called Type 2 
singularities. The famous infrared renormalons of large order perturbation 
theory are connected to singularities of this type. As numerical simulations 
of nonperturbative behavior are necessarily restricted to finite volume lattice 
theories, the question arises whether singularities of this type are really ac- 
cessible using simulation techniques. We shall see below that the question is 
complicated by the nonuniformity of the analytic continuation to positive s 
in the large volume limit, which will require a modification of the definition 



( 15[) in the finite volume situation. For the time being, let us illustrate the 



absence of a true singularity in flTjp in a discrete system with a simple toy 
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model. Consider a mass gap relation connecting a bare coupling / and a 
dynamically generated mass m of the form 

7 = iSsrwi^ ln( — ) ' L ^°° (20) 

The "infinite volume limit" thus gives m = 1 f. , and Borel transforms 
-B(s) of any analytic function G(m) expandable at m = will display simple 
poles at s — 0,1,2,3... For example, if G(m) = jt^, we have (at L = 
oo) B(s) — — 1/s — 1/(1 — s) (for other choices of the "momentum" in the 
denominator of G(m), there will typically be poles at all integer s). On the 
other hand, for finite L, we have 

POO 1 I 

Bl(s) = / e-//— -d(-) 
Jo 1 + m j 

00 1 1 1 L V 1 

(m + n/L) 2 l + m K ' 

where the m-integral in (21) is convergent for Re(s) < 0. We may analytically 
continue to s = 1 by deforming the m contour to avoid a singularity from the 
n = term in the exponent,i.e. by keeping arg(^) fixed near m = 0. This 
can be done without encountering any pinches of the integration contour. 
Consequently, Bl{s = 1) is perfectly finite for finite L. In Section 6 we shall 
see that Bl(s) needs to be modified so that a smooth approach to the infinite 
volume behavior for positive s is obtained. 



3 Non-Borel-summable spin models 
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3.1 Separable models (D=l) 

A simple model which already serves to illustrate several features of the 
renormalon singularity structure in more general non-Borel theories is the 
ID spin chain defined by the action 

S = ^(l-fa-fa+x) (22) 

Here the fa are unit 3-vectors and free boundary conditions are chosen so 
that the model factorizes in a trivial way. Indeed, the partition function is 

Z L {f) = Jli l( ifae- S[;i>]lf 

= ({) L (l-^ 2// ) L (23) 

showing a trivial separation of the "perturbative" (oc f L ) and "nonpertur- 
bative" (oc f L e^ 2n ^) terms. For L — 1, the Borel Transform of Z is 

B 1 (s) = 1 -\n(^) (24) 

with discontinuity 

£i(s) = -Im B(s + ie) = -9{s)6{2 - s) (25) 

7T 2 

Note that poles only appear in B\(s) after a derivative with respect to s. This 
effectively (cf. Eq (|T5|) ) removes the overall perturbative factor of / (associ- 
ated with the measure in the path integral), exposing the nonperturbative 
behavior. For L large, Bl(s) is even smoother (each prefactor of / effec- 
tively integrates with respect to s). The singularity structure is numerically 
invisible in this case. In Fig. ^| we show B 8 (s). It is clearly difficult to tell 
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that there are actually singularities at s = 2, 4, 6, 8, ...! In Section 4 we shall 
see that this perturbative camouflaging of the singularity structure is typi- 
cal of disconnected quantities like Z, but is not a problem when connected 
quantities (with a well-defined infinite volume limit) are considered. 




3 spin chain, 0(3) 
8 plaquette 2D, SU(2) 
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Figure 2: Borel Transform B for L = 8 spin chain 



In agreement with the usual analogy between spin models in D and gauge 
models in 2D dimensions, the separable model of 2D pure lattice gauge theory 
is found to have an essentially identical Borel structure. If we take the gauge 
group to be SU(2), the partition function for a system of L plaquettes (with 
free boundary conditions) is found to be Zi{f) = [2/e _1 /^ii( j)] L . For L — 1, 



the Borel Transform is Bi(s) = J s(2 — s). The Borel discontinuity (=density 
of states in functional integral) for a system of 8 plaquettes is also plotted in 
Fig(^|) for comparison with the spin chain case. Once again, the singularities 
at s = 2, 4, 6, ... are numerically invisible. 
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3.2 Borel Structure of 2D Sigma model 

The partition function of the continuum 2D Sigma model is given by the 
Euclidean functional integral 

Z(f) = J V0(x)e-^f d2x]9 ^ 2 (26) 

Here is a unit N-vector field, / is a bare coupling, and a UV cutoff is 
presumed. On a finite 2 dimensional LxL lattice, the corresponding quantity 
is 

z L (f) = /nA^e^^"^- 1 )-^^ 1 ^ 12 

= (l) NL2 JUnd^dpne^P^NW-^lAM 2 (27) 

Here n labels lattice sites, and an auxiliary condensate field p n has been 
introduced to implement the nonlinear constraint. Integrating out the (f) 
fields 

/ 



Z L (f) = (p NL2 J Tl n dp n e~TSc«{Pn] (28) 



S cS = TrLn(— A 2 - 2ip) + j £ p n (29) 

J n 



where 



For large N, the functional integral in Eq(|28|) is dominated by a translation- 
ally invariant complex saddle point at p n = p sa d where 

7 = 77 £ TTtT > m = -2ip sad (30) 



Here the kinematic momentum factor appearing in the propagator denomi 
nator is d(k) = 4(sin (nk x /L) 2 + sin (nky/L) 2 ). Thus the 1/N expansion [13 
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generates immediately a dynamical (squared) mass m(f). Note that for any 
finite L, the dynamical mass m(f) defined implicitly in Eq (|30|) is actually 
analytic at / = 0. Any analytic function of m (e.g. the large N limit of the 
(j) propagator at momentum k, namely —7^ — ) will therefore have a conven- 
tional Borel transform B(s) which is entire in s. The full Borel transform 
B(s) of the large iV propagator has a cut on the positive axis, but no true 
singularities of B for finite L. The problem of recovering the singularities 
of the large volume limit from lattice data will be discussed in Section 6. 
Type 1 (instanton) singularities which would surface already at finite L are 
of course absent in the large iV limit. 

The large N result for the partition function is 

Z L ~ (f/N) NL *e-^&^i taW ® +m) - m/f) (31) 



Again, as for the spin chain case, the perturbative prefactor will numerically 
camouflage the nonperturbative structure. 

For future reference, we list here several important features of the large N 
limit in the continuum theory. Initially, assume both infrared and ultraviolet 
cutoffs to be present for the modes of the <fi field: A < k 2 < A. The mass gap 
equation then becomes 

1 r d 2 k 1 



/ J (2ir) 2 k 2 + m 



47T A + m 



giving a dynamical squared mass 



A e - A7T/f - A 
m = 1 _ e -W/ (33) 
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displaying the characteristic essential singularity at / = 0. The singularity 
is present even though A 7^ (and conventional perturbative IR renormalons 

— # 

are absent:cf. Section 5)! The large N limit of the propagator in momentum 
space 

1 1 m m 2 

possesses a meromorphic Borel transform B(s). Taking A = for simplicity 
and changing variables to x = 1// 

r+°° Ae~ 4irx 1 n 

/ dxe ax { , - )" = A ra (- + —, + ..) (35) 

Jo v i_ e -47n^ v 4vm-s 47r(n+l)-s ' v 7 

so there are simple poles at s = 4n, 8n, 12ir, .... Switching on the infrared 
cutoff A modifies the residue, but not the location of these poles. This ap- 
parently paradoxical situation is elucidated further in Section 5. 

The meromorphic character of the Borel transform of the propagator 
extends to other Green's functions in the large iV limit. The two-current 
correlator, for example, is defined by 



n""(?) = ^- 2 Jd 2 x<o\Jt; b (x/2)r ab (-x/2)\o> 



J" ab = \<Pady b (36) 

Neglecting terms of 0(q 2 /A.) and subtracting at q — 0, one has (in the large 
N limit) 

n ^ = -^(^-^T) / 1(fem (l + (37) 
o7r q z Jo m 

The Borel transform can be computed by leaving the Feynman parameter 
integration to the end. Again neglecting terms of order q 2 /A (but notm/q 2 \) 
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the relation between m and / is effectively m = Ae 47r ^ so 

B^(s) = Lw 4 -(^_^) f l dx [ A dm m -i-/** In (1 + ^ ~ ^ ) 

327r 2 jo jo m 

(38) 

We are interested in the singularity structure of B^(s) for Re(s) > (the 
generalized Borel Transform is left-half-plane analytic). Since 

1 dx l +0 ° dm m~ l - s '^ In (1 + x(1 ~ ^ ) (39) 
o ja m 

is analytic for Re(s) > 0, we may complete the m integral to J +o ° without 
altering the singularity structure in the right-half-plane. One then finds 

[ ,q) ~32^ [9 J V r(2- S /2vr) [ ' 

which has double poles at s = 47r,87r, ... These poles would be visible as 
sharp peaks if B tJ,v (s) could be computed at Re(s) < and then analytically 
continued to s = sr + 27,7 << sr,sr > 0. Note that B^ u (s) satisfies the 



Brown- Yaffe-Zhai relation [11 

discos; q 2 = -Q 2 ) = sin (nb s)B^(s; Q 2 ), b = — (41) 

47T 

From the foregoing, it is apparent that the nonlinear sigma-model in the 
large N limit exhibits an interesting and elaborate renormalon pole structure 
|10|j in the right-half Borel plane. This is so even though a conventional 
perturbation expansion applied to the large N result leads to trivial results, 
insofar as the nonperturbative structure enters in this limit entirely through 
the quantity m ~ e~ 4n ^ which vanishes to all orders in a formal expansion 
around / = 0. 
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4 Perturbative Camouflage of Renormalon Sin- 
gularities 



Formally the most convenient field-theoretic quantity from the point of view 
of the generalized Borel function B(s) is the basic partition function Z(f) of 
Eq (|13|), as the discontinuity of B(s) in this case is nothing but the density 
of configurations in a given action shell in the functional integral. Unfortu- 
nately, Z does not possess an infinite volume limit, and as we saw in Section 
3, the Type 1 singularities of B(s) for a system of large but finite volume V 
are in any event camouflaged by perturbative prefactors ~ / which smooth 
the s-dependence (as multiplication by / is equivalent to integration by s). 
The Borel transform of Z(f) is very accurately computable by Monte Carlo 
techniques: from the density of states interpretation described in Section 
2 one has only to histogram the frequency of configurations generated in a 
Monte Carlo simulation (relative to a conveniently chosen reference action) 
versus action. The result of a simulation of B(s) for a 8x8 2D sigma model is 
shown in Fig(^), to be compared with Fig(0) for the spin chain or 2D QCD. 
Note that the singularities at positive s are numerically invisible. 

The connected vacuum amplitude W(f) = In Z/V is clearly preferable to 
Z in this regard. Apart from having a well-defined infinite volume limit, per- 
turbative prefactors in Z appear as additive contributions in W, effectively 
uncovering the renormalon structure of the theory. We may expect the same 
to be true of other connected Green's functions (e.g. two-point functions) of 
the theory, which are expressible as source-derivatives of W. For example, 
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consider the trivial spin-chain model of Section (3.1): 

W L = j In Z L = In {f/2) + In (1 - e~ 2 / f ) (42) 

the Borel transform By/ of which exhibits simple poles at s = 2n, n = 1, 2, .. 
with residue 1/n. The continuation to positive s thus exhibits Lorentzian 
peaks: 




32 64 96 128 160 192 224 

S 



Figure 3: Density of configurations for 2D 0(3) Sigma model 

Assuming the Borel discontinuity B(s) of Z has been computed (say, 
by Monte Carlo) sufficiently accurately, the analytic continuation of Z(f) to 
x = j — iy (y positive real) is available as the Fourier transform of B. Recall 

r+oo 

Z(x = 1/f) = / dsB{s)e sx (44) 
Jo 
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while 

r+oo 

B w (s)= dxe sx \nZ(x), Re(s)<0 (45) 
Jo 

Analytically continue B w from real negative s to s = sr + 27 (sr, 7 > 0) by 
simultaneously deforming the x-contour in Eq(^5|) to x — iy. Thus 

lmB w (s R + i-f) = Re / dye lSRy -~ <y In / dsB(s)e lsy (46) 

Jo Jo 




Figure 4: Borel singularities of W in ID 0(3) spin chain 

The result of this computation for the spin chain model is shown in Fig(f|) 
for 7 = 0.05. The singularities at s — 2, 4, 6, .. are clearly visible. (The oscil- 
latory structure visible in Fig([|) is due to the fact that the integral over y, 
performed numerically, was cut off at the upper end). It would be clearly of 
great interest to locate the Type 1 singularities in more interesting lattice field 
theories along these lines. At present most of the discussion of renormalon 
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structure centers on correlation functions (e.g. the Euclidean current-current 
correlator) and we shall focus instead on an alternative approach which al- 
lows the extraction of renormalon singularities (of both types) of connected 
Green's functions. 



5 Status of the First Infrared Renormalon in 



It has been known for some time that evidence for a nontrivial singular- 
ity structure in the right-half-plane of the Borel variable s can be detected 
within the framework of conventional (weak-coupling) perturbation theory. 
Consider an asymptotically free field theory (such as massless QCD) with no 
intrinsic mass scale at the classical level. Any scalar Euclidean momentum 
space correlation function H(Q 2 ) corresponding to a renormalization group 
invariant quantity will have a formal perturbative expansion in powers of the 
running coupling a(Q 2 ): 



An example of great phenomenological interest is the hadronic electromag- 
netic current correlator, whence IT(Q 2 ) = ^ / d 4 xe iqx < 0\T[j^(x)j^(0)]\0 >. 
Perturbative infrared renormalons are associated with subsets of graphs which 
result in a contribution to the coefficients c n at large n of the form 



QCD 



n(Q 2 )~E^« n (Q 2 ) 



(47) 



n 




(48) 



P 
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where 60 — — ilT is the first nontrivial coefficient in the beta-function and 
p is an integer. The behavior ( |4"BD implies a singularity in the naive Borel 
Transform B(s) at s — p/bo. Such terms can be seen to arise from sets 
of graphs generating the leading infrared logarithms in the current-current 
correlator (see Fig(|J). Introducing an infrared cutoff A, the contribution of 
the subset of graphs shown in Figfl5|) is found to be || 

r (y(C) 2 ) /MQ 2 A 2 ) 

/ f d f < ^IA V« n+1 (Q 2 )(M n / x n e~ 2x dx (49) 

h» l + 6 a(Q 2 )ln(J) ^ Jc K J 







^— < 












Q->X 












~— — > 





Figure 5: One of a set of bubble graphs contributing to the leading pertur- 
bative renormalon in QCD 

If no infrared cutoff is present (A = 0), the large n behavior of the x 
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integral is dominated by x ~ n and we have 

c n ~ (|)™n! (50) 

implying a singularity in B(s) at s = The absence of a singularity at 
the first possible renormalon location s = ^ can be related to the absence 
of a local gauge-invariant operator of dimension 2 in QCD (the extra factor 
of I 2 in the numerator of (|49|) is a direct consequence of gauge-invariance). 
The nonperturbative ambiguity associated with the singularity at s = ^ is 
of the form e b o a ^ - 1 ~ i exactly the leading power behavior associated 



with coefficient functions of the lowest dimension nontrivial operator F 



2 

fj,!/ 



appearing in the operator product expansion of this current correlator. 

With an infrared cutoff present, the renormalon singularity disappears, 
as the large n behavior is dominated by x ~ In (^-) and 

Cn~ln(^) (51) 

This power growth means that the contribution from these graphs to B(s) is 
entire in s. Of course, the Taylor series generated by the coefficients (|5T|) has 
a finite radius of convergence (before the Borel transform). The divergence 
occurring when A is decreased to the point where ln(Q 2 /A 2 ) ~ b a(Q 2 ) * s 
just the entrance of the Landau singularity into the range of the momentum 
integration. 

The preceding discussion applies almost word for word to the two current 
correlator ( j36[) in the 2D sigma model. A typical set of diagrams giving 
an infrared renormalon in this theory is shown in Fig(O). Here the current 
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insertions do not decouple at zero momentum and we have, instead of (|49|), 

L dl \+ h ^nZ MM ~ ^ n+1 (Q 2 K r^x^dx (52) 
J\ 2 1 + b a{Q 2 ) In ) J 

If the infrared cutoff is set to zero, the large order behavior is c n ~ b^nl, 
leading to a singularity in B(s) at s — This corresponds to a local operator 
l^^p of dimension 2. If A ^ 0, the perturbative renormalon disappears 
completely, as the large order growth is power rather than factorial. 



Figure 6: One of a set of bubble graphs contributing to the leading pertur- 
bative renormalon in the 2D sigma model 

An important advantage of the 2D sigma model is the availability of the 
large iV saddle-point technique as a reliable tool for exhibiting the mechanism 
of dynamical mass generation in the theory. It was pointed out previously 
that dynamical mass generation occurs even if an infrared cutoff is imposed 
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on the modes of the elementary fields of the theory (recall Q33|), where 
m, A, A represent the squared mass, UV and IR cutoff's resp.). From ( j33|) it 
is apparent that even if A 7^ 0, the Borel Transform B(s) of the propagator 
q2 + m has simple poles at s — Ann, just as in fl35|). However, we have just seen 
that if A 7^ 0, the usual perturbative sources of infrared renormalon singulari- 
ties in the naive transform B(s) disappear! In fact, the renormalon structure 
of the full Borel function B(s) in this case is associated with the appearance 
of a zero-momentum condensate for the auxiliary O(N) singlet field p(x) (see 
(|27|)), which couples to the singlet two <fi channel. Even though the individ- 
ual constituents are restricted to have nonzero momentum k 2 > A, pairs 
of these quanta can condense at zero momentum, precisely as in the case of 
Cooper pairs in BCS superconductivity The moral of this example is clear: 
the full Borel transform B(s) can have singularities of purely nonperturba- 
tive origin, not associated (at least in any direct fashion) with the usual IR 
renormalons of perturbation theory. 

The preceding discussion applies to the renormalon structure of the con- 
tinuum field theory in which the momentum modes are still continuous, but 
a sharp infrared cutoff is introduced. The introduction of a finite lattice im- 
poses an infrared cutoff uniformly on all fields of the theory, including the 
condensate field p, thereby eliminating the singularities at s — Aim of B(s). 

It was pointed out in section 2 that a singularity in B(s) for the Borel 
transform of the partition function Z(f) is directly related to a singularity 
in the density of configurations in the defining functional integral for Z. 
Consider the current-current correlator in quenched QCD (i.e. dropping 
internal quark loops). The UV cutoff will be implemented on a lattice (with 
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a bare coupling a = g 2 /4ir) and the compact link integrations written / T>A^ 
for simplicity. Then (again, representing the lattice Fourier transform in 
continuum notation for simplicity): 

II(Q 2 ) = N(Q 2 ;f)/Z(f) 
N(Q 2 ;f) = ^/d 4 xe^|^Tr[ 7 ^A(x;^) 7 ,A(-x;A M )]e~^] 

= JvA,V[A;q]e-^ (53) 

where A(x; A) is the quark propagator in the background field A. The Borel 
discontinuity of N(Q 2 ; f) is 

Bn(s;Q 2 ) = JvA^P{A ]q ]5{s-S[A^) 

= P(s; q)Bz(s), B z {s) = J VA^s - S[A,}) (54) 

Here V(s; q 2 ) is the average of V[A^; q] over all configurations with action 
equal to s. Bz{s) is the fundamental density of configurations in the theory. 
As pointed out by t'Hooft in his seminal study of the renormalon struc- 
ture of QCD [KJ, Borel singularities in the transform of N(f) or Z(f) will 
generically transfer, via the convolution theorem, to the Borel transform of 
the connected ratio N/Z, barring some miraculous cancellation. In fact, the 
location of renormalon singularities (though not their type and strength) 
are expected to be substantially independent of the specific Green's func- 
tion under consideration: rather, they derive from the basic structure of the 
functional configuration space of the theory. 

It has sometimes been asserted that the presence or absence of renormalon 
singularities can be inferred solely from the assumption of nonperturbative 
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validity of the Wilson operator product expansion. For example, the ab- 
sence of a perturbative IR renormalon in B(s) at s = l/b was related above 
to the absence of local gauge-invariant operators of dimension 2. Such a 
singularity will never be visible in dimensionally renormalized massless per- 
turbative QCD, where the only mass scale in the theory is introduced via 
the regularization procedure in such a way as to result in purely logarithmic 
(rather than power) contributions to renormalized quantities. Rather, we are 
concerned here with intrinsically nonperturbative contributions to the coef- 
ficient function Ci(x) of the identity operator, which is present regardless of 
the subsequent set of nontrivial operators in the theory. Such contributions, 

of the form mi 2 in coordinate space, where m is a dynamically generated 

i_ 

squared mass of the form ~ e b o a , a — > 0, inevitably produce a Borel sin- 
gularity at the first location s — ^. On the other hand, they vanish formally 
to all orders in a weak coupling expansion, and (in a non-Borel theory) are 
not necessarily reconstructible from conventional sources of factorial growth 
in perturbation theory (cf. situation discussed above for 2D sigma model in 
presence of an IR cutoff!). 

Finally, let us note that the presumed existence of a mass gap in the 
theory can be viewed as prima facie evidence for a universal infrared renor- 
malon singularity at the first available location. It is very difficult to see 
how a singularity at s = l/b can be avoided in almost any physical quan- 
tity in QCD once such a mass gap is assumed, together with the renormal- 
ization group and asymptotic freedom. As the simplest example, consider 
the finite temperature partition function of pure QCD at low temperature 
Evidently, for large f3, the one glueball sector dominates and (defining 
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m = (glueball — mass) 2 ) 




l + V 



I 




(2tt) 3 



m = Ae b ° 



(55) 



The higher order terms are suppressed at low temperature, so cannot be 
expected to cancel (except accidentally) any renormalon singularity found 
in the first nontrivial term. The Borel transform f?i(s) of the 1-glueball 
contribution to the free energy InZ/V thus becomes 



A similar singularity is expected to arise in all connected Green's functions 
of the theory once dynamical mass generation occurs. 

Of course, to settle this contentious issue definitively, it will be neces- 
sary to perform reliable, fully nonperturbative calculations of the general- 
ized Borel function B(s). The only tool available for doing this in the case 
of QCD is lattice theory, so we must now face directly the issue of extract- 
ing renormalon structure of field theories formulated on a (necessarily) finite 
space-time lattice. 




(56) 



Integrating by parts to expose the pole at s = -r- we find 




(57) 
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6 Finding renormalon singularities on a finite 
lattice 



If, as we have argued, weak-coupling perturbation theory is not a reliable 
guide to the full renormalon structure in an infrared unstable theory with 
dynamical mass generation, it will be necessary to carry out a fully non- 
perturbative evaluation of the Borel function B(s) for the relevant physical 
quantity. In QCD, this restricts us to a numerical (Monte Carlo) simula- 
tion of the lattice gauge theory. There are two obvious difficulties which 
must be overcome. Firstly, the information obtained in such a calculation 
is necessarily subject to statistical errors. Consequently, any technique em- 
ployed for the extraction of the singularity structure of B(s) for Re(s) > 
(obtained by analytic continuation from Re(s) < 0) must be fairly resistant 
to the inevitably noisy input information. Secondly, simulations can only 
be performed of systems with a finite number of degrees of freedom, i.e. at 



finite volume. The Borel Transform B(s) (15) will in this case lack the Type 
2 singularities associated with the condensation of zero momentum modes 
in the theory. In this section we shall explain one possible modification of 
the definition fll5|) capable of revealing the precursors of the infinite volume 
singularity structure on a finite lattice. 

As usual, a convenient model for studying the renormalon structure on 
a finite lattice is the O(N) sigma model in 2 dimensions: the large N limit 
of the theory is analytically solvable, even on the lattice. We shall use this 
model to gain insight in formulating a modified lattice Borel transform. We 
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begin by pointing out a basic difficulty in the definition (|T5| ) of the full Borel 
function -Bl(s) when the transform is applied to a lattice system with a mass 
gap generated via logarithmic infrared divergencies. To make the algebra as 
transparent as possible, consider the toy model of Section 2 where the gap 
equation is modelled by the discrete sum 

7 = 7E^77^(m)-ln(^), L-oo (58) 
/ L ^ m + n/L m 

Here m is the analog of the squared dynamical mass and we take G(m) = 
p2 \ m as our model of the momentum space propagator, with Borel transform 

B L {s) = / e s "G(m(f))d(-), Re(s) < 
Jo f 

= - / s' L {m)e SSL{m) G{m)dm (59) 
Jo 

This integral defines the analytic function -Bl(s) for Re(s) < 0. The problem 
is that the L —>■ oo limit does not commute with the desired analytic contin- 
uation to the right-half s-plane. This becomes clear if we subject (p9|) to an 
integration by parts 

B L {s) = -U SSl(0) G(0) + -J e SSL{m) G'{m)dm (60) 

For < Re(s) < 1, the integral in ( |60| ) is now well-defined as L — > oo 
(as the integration is convergent for m — > 0), but s_l(0) ~ ln(L) causes a 
power divergence oc L s if the large volume limit is taken after the analytic 
continuation is performed. 

In the case of the sigma model at large N, where we have analytic control 
of the theory, one may devise the following cure for this problem, which arises 
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from the deviation at finite L in the integration measure for small m from 
the infinite volume form, Namely, one devises a modified Borel transform on 
the lattice which has a smooth analytic continuation to Re(s) > at large 
L. In the L —>■ oo limit, the integration measure goes over to m~ 1 ~ s dm. 
Choosing an infrared cutoff (3 ~ 0(1/ L) (the precise value is immaterial) we 
may define an alternative transform of the (nonperturbatively determined!) 
propagator G{m) on the lattice as 

B mod (s) = e vs [ m- l - s G(m - f3)dm + B L (s) (61) 

J 

with T],(3 satisfying r] = sl(0) — ln(l//3), i.e 

e s(s L (0)- V ) = p-s 

The relation (|62|) ensures that the integrand m~ s G(m — f3) in the first term in 



( 6TD matches (at m = /3) smoothly to the corresponding lattice determined 
integrand e SSL ^G(m) in Bl(s). The supplementary contribution generates 
a term — ^e vs G(0) ~ L s as L — > oo, exactly cancelling the pathological term 
in (p0|). We have checked that the above procedure leads to a numerically 
stable detection procedure for renormalon singularities using the test case 
of the 2D sigma model. As input to the Borel transform, the feasibility of 
the method can be studied by using the analytic large N results (instead of 
actual Monte Carlo data). Statistical noise can also be introduced by hand 
to check the robustness of the procedure. The mass gap equation on a finite 
lattice for the 2D sigma model is 

r^?wT^ s i Si(m) (63) 
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where k = {^k x , ^k y ), k x , k y = 1, 2, ..L and d(k) = 4(sin(^) 2 + sin^) 2 )- 
In the large N limit the momentum space <j) propagator is 




Figure 7: Borel Transform of propagator, inverse propagator on a 30x30 
lattice (7 = 0.1, N T = 10) 

In analogy to (|5"T|), we analytically continue the modified transform (sup- 
pressing for notational simplicity the momentum dependence) 

B mod (s) = e" 4 / m- l - s G{m{f))dm - f s' L {m)e s SL ^G(m)dm (65) 
Jo Jo 

with (3 ~ 0(1/L 2 ), r] = sl(Q)— \n(l//3). The analytic continuation to positive 

real s can be performed in a numerically stable way by writing 

s = s R e w (66) 

and Taylor expanding around w = iir (where B mod (s) is analytic). The circle 
of convergence has radius tt, as -B mo d is analytic in the cut plane, and we 
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can evaluate Im B mo a(sR + ij) by summing the series to some finite order 
Nt at w ~ Z7- The result for Nj- = 10,7 = 0.1 on a 30x30 lattice is shown 
in Fig(^) for the Borel transform of both G (with d(q) = 1) and the inverse 
propagator G~ l . The effect of keeping more terms in the Taylor expansion 
is shown in Fig(§) for the Borel Transform of G~ l . The peak at the first 
renormalon becomes sharper as Nt increases. However, the results become 
(for fixed 7) less accurate as sr increases, which requires w to approach the 
edge of the circle of convergence. 




Figure 8: Borel Transform of inverse propagator on a 30x30 lattice (Nt = 
6,10,14) 

Finally, the sensitivity of the continuation procedure to a random su- 
perimposed 1% statistical noise on G(m) (well within the reach of Monte 
Carlo simulations, for example) is shown in Fig(P). The first renormalon is 
still clearly visible at s ~ 1, with the effects of the noise becoming serious 
only for s > 1.5. In all of these cases, we should emphasize that the first 
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term in (|65| ) is absolutely essential to obtain stable results from the analytic 
continuation. 




Figure 9: Borel Transform of noisy propagator on a 30x30 lattice (Nt = 10) 

7 Summary and Discussion 

In this paper we have used the nonlinear sigma model as a convenient ana- 
lytical platform for examining the interrelationship of various versions of the 
Borel transform. Infrared unstable theories (like QCD or the sigma model) 
with dynamical mass generation are inevitably non-Borel summable, and a 
precise Borel reconstruction of such theories is only possible starting from 
the generalized Borel transform B(s), defined as the Laplace transform with 
respect to the inverse bare coupling of the relevant amplitude. (The complex 
nonperturbative behavior of such theories is entirely connected with the in- 
frared structure, so it suffices to work throughout with an ultraviolet cutoff: 
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in the case of QCD, this implies a lattice formulation.) 

Unfortunately, the conventional definition of the Borel transform in terms 
of the formal perturbative series (for which we have used the notation B(s)) 
does not lead to a precise reconstruction theorem in the nonBorel case. This 
function extends (in the sense of having identical asymptotic expansion for 
small s) to a function B(s), the discontinuity of the generalized transform 
B(s) introduced above, which does however yield such a reconstruction. The 
essential difference between the two functions can be seen explicitly in the 
large N limit of the nonlinear sigma model, where an infrared cutoff in the 
elementary modes of the theory can be seen to eliminate the conventional 
infrared renormalons of perturbation theory, while leaving the singularities 
in B(s) intact. 

The necessity for a fully nonperturbative approach to the generalized 
Borel function B(s) leads us to the consideration of simulations on a finite 
lattice, which is often the only tool available for reliable nonperturbative cal- 
culation. The Borel singularities on the positive s axis have been seen to 
fall into two types. Type 1 singularities are present (as actual singularities) 
already in the finite volume theory: however they are numerically camou- 
flaged by perturbative prefactors in disconnected Green's functions. Finding 
such singularities on a finite lattice will therefore require removing discon- 
nected vacuum contributions (cf. Section 4). Type 2 singularities (the usual 
IR renormalons fall into this category) are strictly speaking only present in 
the infinite volume limit. Moreover, the interchange of analytic continua- 
tion to Re(s) > and the infinite volume limit is not uniform, requiring a 
modification (discussed in Section 6) of the definition of the Borel Transform 
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on a finite lattice in order to allow detection of the renormalon precursors. 
Such a modification indeed allows us to extract the location of renormalon 
singularities for the analytically tractable test case of the sigma model at 
large N. The application to situations (say the sigma model with iV = 3, or 
quenched QCD) where Monte Carlo methods must be employed to compute 
the nonperturbative structure of the theory is presently under study. 
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